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SECTION  1 


INTRODUCTION 


1 .1  Simulation  Extensions 

The  objective  of  this  project  is  to  identify  and  quantify  those  generic 
mission,  scenario  and  sensor  parameters  (as  well  as  their  interactions)  which 
drive  the  performance  of  space-based  staring  IR  surveillance  systems.  These 
efforts  have  been  divided  into  two  efforts:  Simulation  Enhancements;  and 
Simulation  Applications. 

Work  in  the  Simulation  Enhancements  area  is  directed  toward  general¬ 
izing  those  models  in  the  Draper  Integrated  Simulations  (DIS)  which  are  too 
limited  in  scope  or  perhaps  restricted  to  one  particular  system.  Models  have 
been  prepared  for  the  platform,  sensor  (including  focal  plane  and  signal  pro¬ 
cessor),  scene,  and  techniques  for  interfacing  the  DIS  with  an  in-house  image 
processor  for  graphics  display. 

The  Simulation  Applications  effort  studies  issues  related  to  generic 
surveillance  system  performance  drivers  using  tools  from  the  Simulation 
Enhancement  effort.  Typical  examples  are  3ensor  line-of-sight  stability, 
focal  D*,  and  critical  scene/scenario  characteristics. 


HALO  Optics 


HALO  is  a  multi-mirror  large  optical  imaging  system  with  a  wide  field 
of  view.  To  maintain  alignment  as  well  as  image  quality,  the  mirrors  are 
actively  controlled  with  actuators.  The  term  "deconvolution"  implies  a  sens¬ 
ing  and  control  scheme  by  which  a  deformed  HALO  optical  system  is  corrected  in 
near  real  time.  In  one  of  the  techniques  worked  on  by  the  Itek  Corporation, 
the  wavefront  errors  of  the  deformed  system  are  measured  with  a  wavefront  sen¬ 
sor  for  several  object  points  near  the  edge  of  the  field-of-view  of  the  sys¬ 
tem.  The  wavefront  sensor  is  placed  around  the  focal  plane.  The  wavefront 
errors  are  decomposed  into  mirror  figure  errors  and  actuator  displacements  are 
generated  to  correct  them  with  a  deconvolution  algorithm.  This  report  de¬ 
scribes  a  limited  blind  test  that  CSDL  prepared  for  Itek  to  test  and  validate 
this  algorithm.  We  conclude  that  the  algorithm  successfuly  determined  the 
actuator  displacements  from  the  wavefront  errors  considered  in  this  test.  The 
residual  wavefront  errors  were  negligibly  small  and  below  the  wavefront  sensor 
noise. 


It  is  recommended  that  one  more  blind  test  be  prepared  for  Itek's 
deconvolution  algorithm  to  include  more  realism.  Some  of  the  desirable  fea¬ 
tures  of  this  test  are  listed  below. 


1 )  Mirror  deformations  will  not  be  produced  by  tho  actuators  used  for 
correction,  because,  in  practice,  a  mirror  is  not  deformed  by  the 
movement  of  actuators. 


2)  Some  static  aberrations  will  be  added  to  the  optical  system  to 
simulate  a  nondiffraction-limited  undeformed  system. 


The  on -axis  wavefront  error  data  may  not  be  excluded,  since  the 
wavefront  sensor  placed  around  the  HALO  focal  plane  cannot  sense 
it. 

Tfte  actuator  influence  functions  used  for  mirror  figure  correction 
will  be  slightly  different  from  those  used  by  Itek  in  their  algo¬ 
rithm,  to  simulate  errors  in  the  knowledge  of  influence  function. 

The  influence  function  for  edge  actuators  may  be  different  from 
that  for  interior  actuators. 


SECTION  2 


SIMULATION  EXTENSIONS  PROJECT 


2.1  Introduction 


The  overall  objective  of  this  project  is  to  identify  and  quantify  those 
generic  mission,  scenario  and  sensor  parameters  (as  well  as  their  interac¬ 
tions)  which  drive  the  performance  of  space-based  staring  IR  surve? llance  sys¬ 
tems.  The  approach  taken  to  achieve  mis  overall  objective  has  proceeded 
along  two  main  avenues:  Simulation  Enhancements;  and  Simulation  Applications. 

Efforts  in  the  Simulation  Enhancements  area  are  directed  toward  general¬ 
izing  those  models  in  the  Draper  Integrated  Simulations  (Dis)  which  are  too 
limited  in  scope  or  perhaps  restricted  to  one  particular  system.  Figure  2-1 
presents  a  simplified  block  diagram  of  the  DIS.  In  support  of  the  Simulation 
Enhancements  effort,  work  has  been  performed  (or  is  in  progress)  c-n  models  for 
the  platform,  sensor  (including  focal  plane  and  signal  processor),  scene,  and 
techniques  for  interfacing  the  DIS  with  an  in-house  image  processor  for 
graphics  display. 

Utilizing  the  tools  developed  as  a  result  of  the  Simulation  Enhancements 
work,  the  Simulation  Applications  effort  studies  issues  related  to  generic 
surveillance  system  performance  drivers.  Typical  examples  are  sensor  line-of- 
sight  stability,  focal  D  ,  and  critical  scene/scenario  characteristics  such 
as  spectral  interval,  orbit,  scene  time-of-day,  etc. 

2 . 2  Simulation  Enhancements 


2.2.1  Platform  Simulation 

Development  of  the  Integrated  Large  Space  Structures  Simulation  (ILS3) 
was  begun  in  1977  in  response  to  a  DARPA  need  for  simulating  the  optical 
performance  of  the  HALO  system.  While  intended  to  be  general  purpose,  the 
internal  structure  of  the  ILS3  (i.e.,  input,  coordinate  systems,  interpola¬ 
tion  methods,  contact  systems,  etc.)  was  rigidly  based  on  the  HALO  design. 

More  recently,  many  analysis  efforts  have  concentrated  on  single  reflector  or 
feedref lector  systems.  In  the  absence  of  a  suitable  raytrace  program  to  pro¬ 
vide  the  necessary  inputs  for  utilizing  the  ILS3,  alternate  special-purpose 
software  packages  were  developed  which  did  not  require  a  raytrace  and  which 
could  compute  wavefront  errors  for  these  simple  paraboloidal  or  spherical  sys¬ 
tems.  While  these  software  packages  were  adequate  for  their  intended  pur¬ 
poses,  this  approach  has  three  major  disadvantages:  the  packages  are  problem- 
specific  and  must  be  modified,  sometimes  extensively,  for  other  applications; 
no  linkage  with  the  DIS  '<?  provided;  and  documentation  tends  to  be  scattered 
and  incomplete. 

The  purpose  of  the  enhancements  to  the  ILS3  is  to  assemble  a  package 
which  can  be  used  to  analyze;  both  multi-mirror  raytraced  systems  and  single 
surface,  paraboloidal  or  eoherical  systems.  This  't  being  done  in  five  steps: 
incorporating  many  of  thj  previously  proposed  ILS  '  changes;  incorporating  a 
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general  purpose  version  of  the  single  mirror  software;  adding  a  compatible 
raytrace  program;  simplifying  data  communication;  and  providing  linkage  to 
the  NAS TP All  and  DISCOS  programs.  The  use  of  NASTRAN  is  important  both  to  pro¬ 
vide  input  to  the  top  level  ILS^  modules  (dynamics,  steady-state,  and  quasi¬ 
static)  and  in  many  cases  to  serve  as  a  replacement  for  them.  These  enhance¬ 
ments  will  provide  significant  new  capabilities  for  dealing  wiJi  both  infrared 
and  radar  systems. 

An  initial  plan  nas  been  formulated  and  preliminary  analysis  completed 
to  accomplish  the  abo'/a  objectives.  The  remainder  of  this  effort  is  cur¬ 
rently  in  process. 

2.2.2  Sensor  Simulation 


Many  aspects  of  the  focal  plane  and  signal  processor  models  originally 
in  the  DIS  were  specific  to  the  Mini-HALO  system.  . In  order  to  be  able  to 
study  effectively  the  capabilities  of  a  generic  space-based  staring  infrared 
surveillance  system,  it  was  necessary  to  generalize  certain  aspects  of  each 
of  these  models. 

The  principal  enhancement  to  the  focal  plane  model  was  to  implement  a 
menu  of  user-selectable  noise  models.  In  addition  to  the  standard  noise  model 
corresponding  to  the  basic  Mini-HALO  focal  plane  (including  saturation  and 
reset  features),  the  user  may  select  that  of  an  ideal  focal  plane  (i.e., 
scene  shot  noise  only)  or  may  specify  either  D*  or  D  /Dg^p  to  characterize 
focal  plane  noise.  In  the  latter  case  the  user  may  further  specify  either 
a  white  noise  or  1/f  noise  spectrum. 

The  initial  signal  processor  model  in  the  DIS  corresponded  to  the  HALO 
Signal  Processor  (HSP).  To  increase  flexibility  and  modularity,  the  software 
architecture  for  the  Signal  Processor  has  been  restructured,  leaving  intact, 
however,  the  core  HSP  algorithm  chain.  In  order  to  be  able  to  deal  with  the 
needs  of  a  general  class  of  space-based  staring  IR  surveillance  systems,  the 
character  of  many  of  these  algorithms  has  been  parameterized,  and  a  number  of 
new  algorithms  have  been  added.  The  focus  of  much  of  this  activity  has  been 
in  the  following  areas;  target  signature  templates;  thresholding  procedures; 
threshold  level  computation  and  scaling;  and  "system  track"  and  "acquired 
track"  criteria. 

The  Sensor  Simulation  Enhancements  are  complete  at  this  point  and  now 
being  used  in  DIS  applications. 

2.2.3  Scene  Simulation 


The  single  most  important  initiative  undertaken  in  the  Simulation 
Extensions  Project  has  been  the  development  of  an  in-house  generic  scene  simu¬ 
lation  in  order  to  bring  within  the  scope  of  the  DIS  the  capability  to  gener¬ 
ate,  manipulate  and  analyze  terrestrial  scene  data  sets  as  functions  of  the 
major  surveillance  system,  scenario  and  mission  parameters.  This  effort  has 
been  accomplished  by  working  in  conjunction  with  Photon  Research  Associates 
(PRA).  Under  subcontract  to  CSDL,  PRA  is  developing  a  software  simulation 
called  GENESSIS  to  meet  DIS  needs  for  scene  generation  and  manipulation.  The 
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PRA  work  is  proceeding  in  two  phases:  GENESSIS-I,  to  provide  a  near-term 
capability  with  certain  simplifying  assumptions  and  a  limited  scene  data  base; 
and  GENESSIS-II,  in  which  all  critical  scene  parameters  will  be  modelled  and 
automated  access  provided  to  the  DMA/Landsat  data  base  to  enable  generation 
of  scenes  corresponding  to  any  region  for  which  DMA  (Defense  Mapping 
Agency)  data  exists#  Tables  2-1  and  2-2  provide  a  summary  of  the  principal 
features  and  current  status  of  GENESSIS-I  and  GENESSIS-II,  respectively. 

Figure  2-2  presents  a  simplified  block  diagram  of  the  GENESSIS  software 
architecture. 

2.3  Simulation  Applications 

The  Simulation  Enhancements described  above  have  considerably  expanded 
the  scope  and  capability  of  the  DIS  *"o  conduct  end-to-end  performance  assess¬ 
ments  of  space-ba?  vd  staring  infrared  surveillance  systems.  As  an  illustra¬ 
tion  of  this  nev  apability,  the  DIS  has  been  used  to  study  the  effect  of 
choice  of  spectral  interval  and  selection  of  scene  local  time  of  day  on  clut¬ 
ter  leakage  through  the  signal  processor.  The  basic  approach  was  to  begin 
with  the  California  coast  (Santa  Cruz)  data  base  and  to  use  GENESSIS  to  pro¬ 
duce  two  scenes  in  the  3. 6- 4.0  pm  band  and  two  scenes  in  the  8. 0-9.0  pm  band, 
each  pair  of  scenes  corresponding  to  8  AM  and  to  noon  local  time,  respec¬ 
tively.  A  nominal  surveillance  system  mission  and  configuration  were  assumed, 
and  the  above  scenes  were  processed  through  the  DIS  under  varying  levels  of 
line-of-sight  (LOS)  jitter.  Tables  2-3  and  2-4  summarize  the  scenario,  sensor 
and  signal  processor  configurations  for  this  investigation.  These  parameter 
values  reflect  the  basic  configuration  for  a  system  whose  principal  objective 
is  the  detection  of  strategic  aircraft  against  a  terrestrial  background.  The 
figure  of  merit  used  to  measure  clutter  leakage  in  the  signal  processor  was 
average  number  of  threshold  exceedances  per  frame.  Figure  2-3  illustrates 
the  four  scenes  used.  The  spectrum  at  the  bottom  of  the  figure  indicates 
relative  temperature  (or  scene  radiance)  within  each  scene. 

Since  the  average  number  of  threshold  exceedances  per  frame  is  a  func¬ 
tion  of  threshold  level  in  the  signal  processor  as  well  as  the  level  of  LOS 
jitter,  the  results  of  this  study  may  be  presented  from  those  two  points  of 
view.  Figure  2-4  shows  the  dependence  on  LOS  jitter  level  of  the  average  num¬ 
ber  of  threshold  exceedances  per  frame.  The  particular  threshold  level 
selected  is  representative  of  what  would  be  used  to  acquire  a  target  with  a 
relatively  high  signal-to-noise  ratio.  The  general  trend  of  the  data  is  such 
that,  for  a  given  level  of  LOS  jitter,  higher  numbers  of  exceedances  are  found 
in  the  LWIR  band  (i.e.,  8. 0-9.0  pm)  than  in  the  SWIR  band  fi.e.,  3. 6-4.0  pm). 
In  either  band,  the  higher  number  of  exceedances  always  occurs  at  noon  local 
time  as  opposed  to  8  AM. 

Figures  2-5  through  2-7  show  the  dependence  of  average  number  of  thresh¬ 
old  exceedances  per  frame  on  threshold  level,  for  progressively  increasing 
values  of  LOS  jitter.  Many  of  the  same  trends  are  also  evident  in  this  data. 
Decreasing  the  threshold  results  in  a  rapid  increase  in  threshold  exceedances. 
At  a  given  threshold  value,  higher  numbers  of  threshold  exceedances  always 
occur  in  the  LWIR  band  as  opposed  to  the  SWIR,  and,  within  each  band  higher 
numbers  of  exceedances  occur  at  noon  in  comparison  with  8  AM.  Increasing  the 
LOS  jitter  level  results  in  a  translation  upward  of  the  curves  at  higher 
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Table  2-1.  Generic  scene  simulation  Phase  I  (GENESSIS-I) 


FEATURES 

•  User-specified  variables 

-Observer  altitude,  zenith  angle 
-Spectral  interval  (2.5-13.0  ym) 

-Atmospheric  model  (LOWTRAN) 

-Fie ld-of- view  location  (within  overall  scene) 

-Scene  angular  resolution 

•  Five  representative  terrain  data  bases 

-California  Coast  (near  Santa  Cruz) 

-Brooks  Range  Mountains  of  Alaska 
-Arctic  Tundra 
-Middle  East 
-Central  Europe 

•  Each  data  base  topology  extracted  from  DMA  data  base 

•  Each  data  base  registered  with  Landsat  data  to  obtain 
materials  assignments 

-Limited  to  14  material  types 

•  Two  representative  cloud  templates  (for  superposition  over 
terrain  data  bases) 

•  Scene  size:  40  km  x  40  km 


STATUS 


Operational  at  CSDL 


Table  2-2.  Generic  scene  simulation  Phase  II  (GENESSIS-II) 


FEATURES 

•  Enhanced  data  base 

-Material  types  expanded  to  20 

-Geometrical  representations  for  three  additional  cloud 
patterns  (Derived  from  NOAA  data) 

•  GENESSIS/DMft'-Landsat  interface 

-Software  package  to  automate  generation  of  input  data  bases 
from  raw  DMA  and  Landsat  data 

-Will  permit  generation  of  any  scenes  for  which  DMA  data 
exists 

•  Generalized  treatment  of  key  atmospheric  effects 

-Continuously  variable  surface-level  parameters  (e.g., 
temperature,  humidity,  wind  velocity,  skyshine,  solar 
scattering) 

•  Optimized  software  architecture 

•  Detailed  user's  manual 


STATUS 

•  GENESSIS-II  (Part  1 )  subcontract  to  PRA  currently  underway 
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Table  2-3.  Scenario  and  sensor  specification  for  GENESSIS 
Santa  Cruz  investigation. 


PARAMETER 

VALUE 

Target 

None 

Backgrounds 

-Location 

Santa  Cruz  (GENESSIS) 

-Spectral  Bands 

3. 6-4.0  pm 

8. 0-9.0  pm 

-Times  of  Day  (Local) 

8  AM 

NOON 

Spacecraft 

-Drift 

None 

-Jitter 

0 . 1  prad  rms 

0.5  prad  rms 

1  .0  prad  rms 

-Altitude 

35,700  km 

-Boresight 

NADIR 

Optics 

-Aperture  Diameter 

1.4  m 

-Obscuration  Ratio 

0.35 

-Thruput 

0.135 

-Point-spread  Function 

Gaussian  (c  =  5.01  prad) 

Focal  Plane 

-Detector  Geometry 

25  x  25  Array  of  100  pm 

(20  prad)  Square  Pixels 

-Frame  Duration 

3.4  sec 

— D* 

2.8  (10)12  cm  /Hz/W 

-Responsivity  non-uniformity 

.Fixed  Pattern 

0% 

.Random 

15% 

0 


Table  2-4.  Signal  processor  configuration  for  GENESSIS 
Santa  Cruz  investigation. 


PARAMETER 

VALUE 

Difference  Filter 

3rd  Order 

Threshold  Level 

User  Input 

Assumed  Target  Contrast 

Negative 

Threshold  Templates 

CSDL  MTD 

Threshold  Algorithm 

CSDL  Rationalized 
Exceedance  Thresh¬ 
olding  (Version  1) 

Cluster  Basis 

Threshold  Exceedances 

Centroid  Basis 

Threshold  Exceedances 

System  Track  Criterion 

3  out  of  3 

Acquired  Track  Criterion 

10  out  of  12 

System  Track 

Threshold  Scaling 

0.5 

Least  Squares  Track 
Linearity  Criterion 

±1  .8  Pixels 

AVERAGE  NUMBER  OF  EXCEEDANCES  PER  FRAME 


Figure  2-4.  Average  clutter  leakage  as  a  function  of  line-of-sight 
jitter  level  for  the  GENE3SIS  Santa  Cruz  scene,  with 
local  time-of-day  and  spectral  band  as  parameters. 


Figure  2-5,  Average  clutter  leakage  as  a  function  of  threshold  level 
for  the  GSrJESSIS  Santa  Cruz  scene,  with  local  time  of  day 
and  spectral  band  as  parameters.  (LOS  Jitter  =0.1  prad 


iH 


|  X  -  8.0  -  9.0  pm 


LOS  JITTER 
FRAME  SIZE 

NUMBER  OF  DATA  FRAMES 
FRAME  DURATION 


=  0.5  prud  rms 
=  25  X  25  PIXELS 
=  25 
=  3.4  s 


X  =  3.6  -  4.0  pm 


THRESHOLD  VALUE 


Figure  2-6.  Average  clutter  leakage  as  a  function  of  threshold  level 
for  the  GENESSIS  Santa  Cruz  scene  with  local  time  of  day 
and  spectral  band  as  parameters.  (LOS  Jitter  =0.5  prad 
RMS) 


AVERAGE  NUMBER  OF 
EXCEEDANCES  PER  FRAME 


X  ■  8.0  -  9.0 


THRESHOLD  VALUE 


Figure  2-7.  Average  clutter  leakage  as  a  function  of  threshold  level 
for  the  GENESSIS  Santa  Cruz  scene,  with  local  time  of  day 
and  spectral  band  as  parameters.  (LOS  Jitter  =  1.0  lirad 
RMS) 
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threshold  values,  as  can  be  seen  by  comparing  Figures  2-5  and  2-6*  In  addi¬ 
tion,  there  is  a  compression  within  each  spectral  band  of  the  exceedance 
curves  corresponding  to  noon  and  to  8  AM*  In  Figure  2-7,  an  LOS  jitter  level 
of  1.0  yrad  rms  (ten  times  the  value  used  in  Figure  2-5)  results  in  the  co¬ 
alescing  of  these  two  curves  for  the  LWIR  band. 

A  subject  of  great  interest  in  the  study  of  signal  processor  performance 
and  scene/sensor  interaction  is  a  determination  of  the  conditions  under  which 
false  target  tracks  arise.  Figure  2-8  summarizes  graphically  the  false  tracks 
which  resulted  during  the  course  of  the  present  study.  Separate  numbers, 
i.e.,  "1"  versus  "2”,  are  used  to  designate  the  relative  position  of  distinct 
tracks  within  the  overall  25  x  25  pixel  field-of-view.  Consistent  with  the 
data  presented  in  Figures  2-4  thi  »ugh  2-7,  false  tracks  were  found  at  low 
threshold  levels  where  the  average  number  of  exceedances  is  very  high.  Under 
these  conditions  the  likelihood  of  a  time  sequence  of  individual  exceedances 
occurring  in  such  a  fashion  as  to  satisfy  the  signal  processor  track  formation 
and  acquisition  criteria  is  considerably  increased. 

2 . 4  Conclusions 

The  focal  plane  and  signal  processor  simulation  enhancements  are  com¬ 
plete,  and  the  platform  simulation  enhancements  are  in  process  and  proceeding 
satisfactorily.  Phase  I  of  the  Generic  Scene  Simulation  (GENESSIS)  activity 
is  complete  and  operational  at  CSDL.  The  last  phase  of  this  activity  is  cur¬ 
rently  in  process  under  subcontract  to  CSDL  at  Photon  Research  Associates. 

The  Draper  Integrated  Simulations  thus  represent  a  highly  sophisticated  tool 
for  end-to-end  performance  trade-offs  of  a  variety  of  space-based  systems,  and 
especially  staring  mosaic  IR  surveillance  systems. 

LOS-induced  clutter  leakage  through  the  signal  processor  has  been  simu¬ 
lated  as  a  function  of  spectral  interval  and  scene  time-of-day  for  a  generic 
staring  mosaic  IR  surveillance  system.  The  functional  dependence  of  the 
scene/sensor  interaction  is  complex  and  has  been  described  quantitatively  by  a 
series  of  curves  showing  average  number  of  threshold  exceedances  per  frame 
both  as  a  function  of  LOS  jitter  for  a  particular  threshold  and  as  a  function 
of  threshold  for  particular  LOS  jitter  values.  Qualitatively  it  is  clear  that 
increased  LOS  jitter  yields  increased  clutter  leakage  in  all  cases  consid¬ 
ered.  Clutter  leakage  in  the  short-wave  infrared  (SWIR)  band  of  3. 6-4.0  ym  is 

significantly  less  than  in  the  long-wave  infrared  (LWIR)  band  of  8. 0-9.0  ym. 

For  either  SWIR  or  LWIR,  clutter  leakage  at  noon  (local  time)  is  greater  than 

early  morning.  For  a  fixed  threshold  value,  increased  LOS  jitter  drives  early 
morning  clutter  levels  toward  noon  clutter  levels  in  both  the  SWIR  and  the 
LWIR  bands.  And,  finally,  low  threshold  values  in  the  presence  of  LOS  jitter 
can  result  in  the  generation  of  false  acquired  tracks. 
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SECTION  3 


HALO  OPTICS 
DECONVOLUTION  TEST 


3 . 1  Introduction 


HALO  is  a  multi-mirror  large  optical  imaging  system  with  a  wide  field 
of  view.  To  maintain  alignment  as  well  ac.  image  quality,  the  mirrors  are 
actively  controlled  with  actuators.  The  teim  "dec...: volution"  implies  a 
sensing  and  control  scheme  by  which  a  deformed  HALO  optical  system  i-~  cor¬ 
rected  in  near  real  time.  In  one  of  the  techn  ques  worked  on  by  the  Itek 
Corporation,  the  wavefront  errors  of  the  deformed  system  are  measured  with  a 
Wcvefront  sensor  for  several  object  points  near  the  edge  of  the  field-of- 
view  of  the  system.  The  wavefront  sensor  is  placed  around  the  focal  plaae. 
The  wavefront  errors  are  decomposed  into  mirror  figure  errors  and  actuator 
displacements  are  generated  to  correct  them  with  a  deconvolution  algorithm. 
This  report  describes  a  limited  blind  test  that  CSDL  prepared  for  itek  to 
test  and  validate  this  algorithm.  We  conclude  that  the  algorithm  success¬ 
fully  determined  the  actuator  displacements  from  the  wavefront  errors. 

3 . 2  Optical  System  for  Deconvolution  Test 

For  a  hardware  demonstration  of  their  deconvolution  algorithm,  Itek 
has  been  planning  to  use  an  optical  imaging  system  consisting  of  three 
active  flat  mirrors  and  a  pair  of  lens  doublets  mounted  back-to-back  in  a 
1:1  magnification.  A  schematic  of  this  system  is  shown  in  Figure  3-1. 

The  lenses  are  the  imaging  elements  of  the  system,  and  the  mirrors 
are  the  deformable  elements.  Aberrations  are  introduced  into  the  system  by 
deforming  the  mirrors.  The  active  mirrors,  thus  simulate  the  deformable 
aspects  of  a  HALO  system.  The  footprints  on  the  mirrors  of  rays  from  a 
given  point  object  are  similar  to  those  on  the  S-HALO  or  WALRUS  mirrors.  A 
perspective  view  o.f  the  optical  system  is  given  in  Figure  3-2.  It  shows  the 
coordinate  system  chosen  to  designate  the  actuator  and  ray  distributions. 

3.3  Actuator  Distribution  and  Influence  Function 


Each  of  the  three  mirrors  has  37  actuators  distributed  on  an  equi¬ 
lateral  triangular  lattice  as  shown  in  Figure  3-3.  In  consultation  with 
Itek,  the  influence  function  of  the  actuators  was  chosen  to  be  similar  to 
the  one  for  active  mirrors  developed  for  the  compensated  imaging  program. 
This  influence  function  is  given  by 


f (r )  =  [l  -  1 .2(r/d) 2  ] 


r  <  d//i72 


»  0  ,  r  >  d/ZTTT  , 


/ 


where 
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r  =  distance  from  the  actuator, 
d  =  actuator  spacing. 

The  variation  of  f(r)  with  r  is  shown  in  Figure  3>-4.  It  is  noted 
that  the  influence  function  is  radially  symmetrical  about  the  actuator 
position.  Moreover,  one  actuator  does  not  produce  any  influence  on 
another.  The  factor  of  1.2  was  chosen  to  minimize  the  effect  of  quilting 
when  all  actuators  are  displaced  by  the  same  amount.  With  a  1.2  factor, 
quilting  is  at  most  8%.  When  an  actuator  is  displaced  by  an  amount  h,  the 
deformation  produced  in  the  mirror  is  given  by  hf(r).  In  our  problem  d  = 
3cm. 


3.4  Mirror  Deformations 


The  deformations  of  the  mirrors  used  in  the  deconvolution  test  are 
given  in  Figure  3-5.  The  numbers  in  this  figure  give  the  actuator 
displacements  in  units  of  X/100.  Mirror  M^  was  deformed  according  to 
cryogenic  deformations  of  an  Itek  HALO  mirror.  Mirror  M2  has  deformations 
described  by  Zernike  polynomials  representing  coma  and  astigmatism.  Mirror 
M3  was  given  a  rigid-body  displacement  along  its  normal.  The  mean, 
peak-to-peak,  and  standard  deviation  of  the  actuator  displacements  are 
indicated  in  the  figure  for  each  mirror. 

3.5  Selection  of  Point  Objects 

Nine  point  objects  were  selected  for  ray  tracing  and  determining  the 
system  wavefront  errors.  Eight  of  these  are  uniformly  distributed  around 
the  edge  of  the  field-of-view  of  the  system  (1.5°  from  the  axis),  and  one  is 
placed  on  the  system  axis  as  indicated  in  Figure  3-6.  The  on-axis  point 
object  was  used  only  for  the  sake  of  completeness  to  include  the  deforma¬ 
tions  of  the  central  portion  of  mirror  Mj.  The  wavefront  sensors  placed 
around  the  focal  plane  cannot  sense  the  effect  of  these  deformations.  In 
practice,  the  central  portion  of  M^  will,  for  example,  be  corrected  by 
CSDL's  image  sharpening  approach.  ' 

3.6  Mirror  Footprints 

Figure  3-7  shows  *'he  footprints  on  the  mirrors  of  the  image  forming 
rays  from  point  objects  1  and  2.  We  note  that,  as  in  a  HALO  optical  system, 
the  footprints  of  rays  from  point  objects  such  as  2,  near  the  edge  of  the 
field  of  view,  on  mirrors  Mj  and  M£  cover  the  mirror  portions  used  in  imag¬ 
ing.  However,  their  footprint  on  mirror  Mj  is  not  only  small  compared  to 
the  mirror  size,  but  also  it  is  off  the  mirror  center.  Accordingly, 
the  footprints  of  the  rays  from  the  8  off-axis  point  objects  cover  the  outer 
portion  of  the  mirror  M}  but  not  its  central  portion. 

The  image  of  the  on-axis  point  object  is  formed  at  the  center  of  the 
focal  plane.  The  central  portion  of  mirror  M^  can  be  corrected  by  the  image 
sharpening  technique  using  the  focal  plane  data. 


Figure  3-5b.  Mirror 
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3.7 


Wavefront  Errors 


When  the  optical  system  deforms,  the  wavefront  error  introduced  at  a 
certain  point  on  the  wavefront  is  given  by  the  change  in'  the  path  length  of 
a  ray  passing  through  that  point.  The  change  in  the  path  length  of  a  ray 
from  a  given  point  object  in  terms  of  the  mirror  deformations  is  given 
according  to3 


AWj.  = 


3 

I 

m=1 


n. 

lm 


Sr. 

im 


cos  0 


im' 


where 

A 

n.  =  unit  vector  along  the  normal  to  the  surface  of  undeformed  mirror 
m  at  the  point  of  incidence  of  the  ray  i, 

5r^m  3  displacement  of  the  point  of  incidence  of  ray  i  on  mirror  m, 

0^„n  3  angle  of  incide  ice  of  ray  i  on  the  undeformed  mirror  m. 

Since  we  are  dealing  with  flat  mirrors,  nim  for  a  given  mirror  is  indepen¬ 
dent  of  i. 

Figure  3-8  shows  the  distribution  of  89  rays  in  the  pupil  used  for 
determining  the  wavefront  errors.  The  wavefront  errors  corresponding  to 
point  objects  1  and  2  when  the  optical  system  is  deformed  according  to 
mirror  deformations  of  Figure  3-5  are  given  in  Figure  3-9.  These  wavefront 
errors  include  X/20  random  error  to  simulate  the  noise  of  a  wavefront 
sensor.  We  provided  Itek  .with  noisy  wavefront  errors  at  89  points  for  each 
of  the  nine  point  objects. 

3.8  Test  Results 


3.8.1  Mirror  Deformations 


From  the  wavefront  error  data  for  the  nine  point  objects  we  provided 
to  Itek,  they  determined  the  actuator  corrections  (displacements)  required 
:o  correct  the  figure  of  the  mirrors  to  reduce  the  wavefront  errors.  The 
actuator  corrections  determined  by  their  deconvo.lution  algorithm  are  given  . 
in  Figure  3-10.  The  residual  errors  obtained  by  subtracting  these  correc¬ 
tions  from  the  errors  given  in  Figure  3-5  are  also  shown  in  Figure  3-10. 

The  mean  and  the  standard  deviation  of  the  residual  displacements  are  listed 
in  the  figure  for  each  of  the  three  mirrors.  We  note  that  the  residual 
errors  are  quite  small  except  for  the  corner  actuators.  The  region  of 
mirrors  where  corner  actuators  are  located  is  not  used  for  imaging.  Hence, 
the  effect  of  these  actuators  on  the  mirror  figure  is  not  observed  by  the 
wavefront  sensor. 
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3.8.2  Wavefront  Errors 


The  wavefront  errors  of  the  corrected  system  were  determined  by 
tracing  89  rays  as  in  the  case  of  the  deformed  system.  These  errors  were 
quite  small  (c^  <  X/40)  indicating  an  excellent  correction.  Random  noise 
of  X/20  was  added  to  these  errors  to  simulate  the  wavefront  errors  that  a 
wavefront  sensor. would  measure.  The  noisy  residual  wavefront  errors  for 
point  objects  1  and  2  are  shown  in  Figure  3-11.  The  mean,  peak-t.o-peak ,  and 
standard  deviation  of  the  wavefront  errors  are  also  given.  It  is  evident 
from  the  standard  deviation  that  most  of  the  error  is  indeed  due  to  the 
wavefront  sensor. 

3.8.3  Point-Spread  Functions 


The  aberrated  and  the  corrected  point-spread  functions  for  point 
objects  1  and  2  are  shown  in  Figures  3-12  and  3-13,  respectively.  The 
standard  deviation  of  the  aberration,  the  Strehl  ratio  (aberrated  irradiance 
at  the  center  relative  to  its  unaberrated  value)  and  peak  irradiance  are 
also  given  in  the  figures.  For  example,  for  point  object  1,  the  deformed 
system  gives  a  Strehl  ratio  of  only  0.025.  The  light  is  spread  over  a  region 
31 XF  wide  compared  with  approximately  2XF  for  the  corrected  (or  aberration- 
free)  system.  With  a  correction  determined  by  Itek  deconvolution  algorithm, 
a  Strehl  ratio  of  approximately  0.95  is  obtained.  The  corrected  Strehl 
ratios  in  Figures  3-12  and  3-13  are  slightly  lower  because  of  the  wavefront 
sensor  noise. 

3.9  Discussion,  Conclusions,  and  Recommendations 


An  optical  system  using  three  flat  mirrors  as  deformable  elements  has 
been  used  for  testing  the  deconvolution  algorithm  developed  by  Itek.  Each 
mirror  has  37  actuators,  and  89  rays  have  been  used  to  determine  wavefront 
errors  of  a  deformed  system.  CSDL  provided  Itek  with  a  known-deformation 
test  and  a  blind  test.  The  known-deformation  test  was  prepared  so  that  CSDL 
and  Itek  softwares  could  be  tested  for  compatibility.  On  the  blind  test, 
two  of  the  9  wavefronts  that  Itek  was  given  are  shown  in  Figure  3-9. 

Itsk's  deconvolution  algorithm  determined  the  actuator  displacements 
on  each  mirror  and  these  were  subtracted  from  the  input  displacements  given 
in  Figure  3-5.  The  corrected  system  was  ray-traced  to  determine  its  imaging 
quality.  It  gave  a  Strehl  ratio  of  approximately  0.9  indicating  an 
excellent  correction. 

We  conclude  that  Itek's  deconvolution  algorithm  successfully 
determined  the  actuator  displacements  from  the  wavefront  errors.  The 
residual  wavefront  errors  were  negligibly  small  and  below  the  wavefront 
sensor  noise. 

We  recommend  that  one  more  blind  test  be  prepared  for  Itek's  decon¬ 
volution  algorithm  to  include  more  realism.  Some  of  the  desirable  features 
of  this  test  are  listed  below. 


are  the  optical 


1 .  Mirror  deformations  will  not  be  produced  by  the  actuators  used 
for  correction,  because,  in  practice,  a  mirror  is  not  deformed  by 
the  movement  of  actuators. 

•  .  % 

2.  Some  static  aberrations  will  be  added  to  the  optical  system  to 
simulate  a  nondiffraction-limited  undeformed  system. 

3.  The  on-axis  wavefront  error  data  may  not  be  provided,  since  the 
wavefront  sensor  placed  around  the  HALO  focal  plane  can  not  sense 
it. 

4.  The  actuator  influence  functions  used  for  mirror  figure 
correction  will  be  slightly  different  from  those  used  by  Itek  in 
their  algorithm,  to  simulate  errors  in  the  knowledge  of  influence 
function . 

5.  The  influence  function  for  edge  actuators  may  be  different  from 
that  for  interior  actuators. 

6.  A  segmented  mirror  will  be  used. 

7.  The  mirror  size  will  be  approximately  the  same  as  the  illuminated 
region . 
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«  MISSION  5 

1  of  I 

§  •  Rome  Air  Development  Center  2 

^  RAPC  p£ana  and  executes  toeatcli,  development,  test  and  ^ 

'  detected  ac.qiuAiti.on  pfiogfiams  in  support  oi  Command,  Control  Q 

2  Communications  and  Intelligence  (C3 7)  activities.  Technical  ? 

i  englnzertng  support  mthin  oJieaA  oi  technical  competence  Q 
S  AA  provided  to  CSV  Ptogaam  OHices  (POs)  and  otkeA  ESV  § 

^  elements.  T^e  principal  technical  mission  aaeas  atie  \ 

.  communications ,  electromagnetic  guidance  and  control,  su/l-  a 
J  yet££anee  oi  ground  and  aerospace  objects,  intelligence  data  i 
»  collection  and  handling,  incarnation  system  technology,  Q 

s  ^ono^pn^u.c  pfiopaga£conp  bo&jL  6toutz  6c, icnceA,  itkoaouxlvz 
t  phyt^QA  and  eZzct/ioyicc  KoZLabltity,  wcUntainab^ititu  and 
S  compatibility. 


